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Abstract 

We study the detectability of the lightest CP-odd Higgs boson of the NMSSM, ai, 
at the LHC through its production in association with a bottom-quark pair followed 
by the a\ — > bb decay. It is shown that, for large tan /3 and very high luminosity of 
the LHC, there exist regions of the NMSSM parameter space that can be exploited 
to detect the a\ through this channel. This signature is a characteristic feature of 
the NMSSM in comparison to the MSSM, as a\ masses involved are well below those 
allowed in the MSSM for the corresponding CP-odd Higgs state. 

1 Introduction 

Supersymmetry (SUSY) is one of the preferred candidates for physics beyond the Standard 
Model (SM). The simplest version of SUSY is the Minimal Supersymmetric Standard Model 
(MSSM). This realisation however suffers from two critical flaws: the /i-problem and the 
little hierarchy problem. The former one results from the fact that the superpotential has a 
dimensional parameter, /x (the so-called 'Higgs(ino) mass parameter'), that, due to SUSY, 
its natural values would be either or the Plank mass scale; yet, phenomeno logically, in 
order to achieve Electro- Weak Symmetry Breaking (EWSB), it is required to take values 
of order of the EW scale, 100 GeV, or, possibly, up to the TeV range. The latter one 
emerges from LEP, which failed to detect a light CP-even Higgs boson, thereby imposing 
severe constraints on a SM-like Higgs boson mass, 114 GeV being its lower limit from data, 
thereby requiring unnaturally large higher order corrections from both the SM and SUSY 
particle spectrum (chiefly, from third generation quarks and squarks) in order to pass such 
experimental constraints (recall that at tree level the lightest CP-even Higgs boson mass of 
the MSSM is less than M z ). 
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The simplest SUSY realisation beyond the MSSM which can solve these two problems 
at once is the Next-to-Minimal Supersymmetric Standard Model (NMSSM) pp. This model 
includes a singlet superfield in addition to the two MSSM-type Higgs doublets, giving rise to 
seven Higgs bosons: three CP-even Higgses ^1,2,3 ( m h! < m h 2 < m h 3 ), two CP-odd Higgses 
01,2 ( m ai < m a 2 ) an d a P a i r °f charged Higgses h^. When the scalar component of the singlet 
superfield acquires a Vacuum Expectation Value (VEV), an 'effective' /i-term, fi e g, will be 
automatically generated and can rather naturally have values of order of the EW/TeV scale, 
as required. Moreover, the NMSSM can solve the little hierarchy problem as well, or at least 
alleviate it greatly, in two ways: firstly, a SM-like Higgs boson can unconventionally decay 
into two ai's with m ai < 2rrib [2], thus avoiding current Higgs bounds (yet this mass region 
is highly constrained by ALEPH [3] and BaBar j3] data); secondly, a CP-even Higgs (h\ or 
/12) has naturally reduced couplings to the Z boson due to the additional Higgs singlet field 
of the NMSSM and the ensuing mixing with the Higgs doublets. 

One of the primary goals of present and future colliders is looking for Higgs bosons. In 
regard to the Higgs sector of the NMSSM, there has been some work dedicated to explore 
the detectability of at least one Higgs boson at the Large Hadron Collider (LHC) and the 
Tevatron. In particular, some efforts have been made to extend the so-called 'no-lose theorem' 
of the MSSM - stating that at least one Higgs boson of the MSSM should be found at the 
LHC via the usual SM-like production and decay channels throughout the entire MSSM 
parameter space [5] - to the case of the NMSSM [S1EIIH]- By assuming that Higgs-to-Higgs 
decays are not allowed, it was realised that at least one Higgs boson of the NMSSM will be 
discovered at the LHC. However, this theorem could be violated if Higgs-to-SUSY particle 
decays are kinematically allowed (e.g., into neutralino pairs, yielding invisible Higgs signals) 

EHUD]. 

So far, there is no conclusive evidence that the 'no-lose theorem' can be confirmed in the 
context of the NMSSM. In order to establish the theorem for the NMSSM, Higgs-to-Higgs 
decay should be taken into account, in particular hi — > a\a\. Such a decay can in fact be 
dominant in large regions of NMSSM parameter space, for instance, for small Ak [H], and 
may not give Higgs signals with sufficient significance at the LHC. 

Besides, there have also been some attempts to distinguish the NMSSM Higgs sector from 
the MSSM one, by affirming a so-called 'more-to-gain theorem' [T2J [THJ [HJ [TTJ [TJ)]. That is, 
to assess whether there exist some areas of the NMSSM parameter space where more and/or 
different Higgs bosons can be discovered at the LHC compared with what is expected from 
the MSSM. Some comparisons between NMSSM and MSSM phenomenology, specifically in 
the Higgs sectors of the two SUSY realisations, can be found in [16J. 

In this analysis, we explore the two theorems at once through studying the direct pro- 
duction of a very light a\ (with m ai < Mz) in association with 6-quark pairs followed by the 
gji — > bb decay, hence a 46-quark final state, at the LHC. This channel has a very large cross 
section at large tan/3 yet, being a totally hadronic signal, is plagued by very large (both 
reducible and irreducible) backgrounds. 

This work is complementary to the one carried in [TTJ [T3], in which we explored the 
t + t - , 77 and decay modes of such a light cii state (again, produced in association 

with 6-quark pairs). 
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This paper is organised as follows: in Sec. 2, we describe the parameter space scan per- 
formed and give inclusive event rates for the signal. In Sec. 3, we analyse signal and QCD 
backgrounds for some benchmark points. Finally, we summarise and conclude in Sect. 4. 



2 Parameter Scan and Inclusive Signal Rates 

In our exploration of the Higgs sector of the NMSSM, we used here the fortran package 
NMSSMTools developed in Refs. [171 [18]. This package computes the masses, couplings and 
decay rates (widths and Branching Ratios (BRs)) of all the Higgs bosons of the NMSSM in 
terms of its parameters at the EW scale. NMSSMTools also takes into account theoretical as 
well as experimental constraints from negative Higgs searches at LEP [TJTj and the Tevatron3, 
along with the unconventional channels relevant for the NMSSM. 

The features of the scan performed have been already discussed in Refs. [TTJ [15], to 
which we refer the reader for details. We map the NMSSM parameter space in terms of six 
independent input quantities: the Yukawa couplings A and k, the soft trilinear terms A\ and 
A K , plus tan/3 (the ratio of the VEVs of the two Higgs doublets) and /i e ff = X(S) (where (S) 
is the VEV of the Higgs singlet). 

For successful data points generated in the scan, i.e., those that pass both theoretical and 
experimental constraints, we used CalcHEP [20] to determine the cross-sections for NMSSM 
Higgs production^. As the SUSY mass scales have been arbitrarily set well above the EW one 
(see Refs. [Tij [TB]). the NMSSM Higgs production modes exploitable in simulations at the 
LHC are those involving couplings to heavy ordinary matter only. Amongst the production 
channels onset by the latter, we focus here on gg,qq —> bb a±, i.e., Higgs production in 
association with a 6-quark pair. This production mode is the dominant one at large tan/3. 
We chose gg,qq — > bbai — > bbbb also because 6-tagging can be exploited to trigger on the 
signal and enable us to require four displaced vertices in order to reject light jets. The ensuing 
46 signature (in which we do not enforce a charge measurement) has already been exploited 
to detect neutral Higgs bosons of the MSSM at the LHC and proved useful, provided tan (3 
is large and the collider has good efficiency and purity in tagging 6-quark jets, albeit for the 
case of rather heavy Higgs states (with masses beyond Mz, typically) (22], [23] . 

As an initial step of the analysis, we have computed the fully inclusive signal production 
cross-section times the decay BR against each of the six parameters of the NMSSM. Figs. 
1 and 2 present the results of our scan, the first series of plots (in Fig. 1) illustrating the 
distribution of event rates over the six NMSSM parameters plus as a function of the BR 
and of m ai . The plots in Fig. 2 display instead the correlations between the a± — > bb decay 
rate versus the a x mass and the ai — > 77 decay rate. It is clear from Fig. 1 that, for our 
parameter space, the large tan/3 and small /i e ff (and, to some extent, also small A) region 
is the one most compatible with current theoretical and experimental constraints, while the 
distributions in n, A\ and A K are rather uniform (top six panes in Fig. 1). From a close look 

1 Speculations of an excess at LEP which could be attributed to NMSSM Higgs bosons are found in [2]. 
2 We adopt herein CTEQ6L [21] as parton distribution functions, with scale Q — the centre-of-mass 
energy at parton level, for all processes computed. 
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at the bottom- left pane of Fig. 1, it is further clear that the BR(ai — > bb) is dominant for 
most points in the parameter space, about 90% and above. In addition, by looking at the 
the bottom-right pane of Fig. 1, it is remarkable to notice that the event rates are sizable 
in most regions of parameter space, topping the 10 7 fb level for small values of m ai and 
are decreasing rapidly with increasing m ai . However, there are some points in parameter 
space, with m ai between 40 to 120 GeV, as shown in the left pane of Fig. 2, in which the 
BR(ai — > bb) is reduced due to the enhancement of the BR(ai — >■ 77) (see right pane of the 
same figure), phenomenon peculiar to the NMSSM and which depends upon the amount of 
Higgs singlet- doublet mixing, see [TTj . 

3 Signal-to-Background Analysis 

We perform here a partonic signal-to-background (S/B) analysis, based on CalcHEP results. 
We assume y/s = 14 TeV throughout for the LHC energy. We apply the following cuts in 
our calculations: 

AR(j,j)>0A, \ V (j)\<2.5, P T (j)> 15 GeV. (1) 

Here, we assume that the 6-tagging probability of a 6-quark is 50% while the mis-tagging 
probability of a gluon is 1% and the one for a light quark is 2%. Figs. 3-7 show the distri- 
butions of the invariant mass of a two 6-jet (di-jet) system after multiplying the production 
times decay rates (after cuts) by the aforementioned efficiency/rejection factors, i.e., true 
6-tagging and mis-tagging probability. It is clear that the largest background is the irre- 
ducible one bbbb, which is one order of magnitude larger than the reducible background bbgg. 
Further, the other reducible background, involving light quarks, labelled here as bbcc, is neg- 
ligible compared to the other two. Notice that all these backgrounds reach their maximum 
in invariant mass at around 40 GeV. Our plots have a bin width of 1 GeV and account for 
all combinatorial effects (as appropriate in absence of jet-charge determination). 

To claim discovery of the lightest CP-odd Higgs, ai, at the LHC, we plotted in Fig. 8 
both signal significances, left-panes, and corresponding signal event rates, right-panes, as a 
function of the collected luminosity, and after integrating over 10 and 5 bins (thus mimicking 
a more optimistic and a more conservative, respectively, detector resolution in mass of a di- 
jet system). From this figure, there is a small window to detect a\ with masses between 35 
at 50 GeV at a very large luminosity, L = 300 fb" 1 , which could well pertain to the final 
number of recorded events of the LHC at design luminosity. Further, for an upgraded LHC, 
with a tenfold increase in design luminosity, known as the Super-LHC (SLHC), which could 
well lead to data sample with L = 1000 fb~ x , there is clear potential to detect an a\ with 
mass up to 80 GeV or so. This is the more probable the better the experimental resolution, 
naturally. 

Finally, notice in Figs. 9-10 the distributions in pseudorapidity-azimuth (the standard 
cone measure) separation of di-jet pairs and average transverse momentum of the jets, re- 
spectively (in the latter case we also distinguish, in the reducible backgrounds, between 6-jets 
and non-6-jets, by exploiting our knowledge of the Monte Carlo 'truth'). Guided by Refs. 
[22| [23] . this has been done to check whether further cuts (in addition to those in eq. ([!])) 
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could be employed to improve the chances of detecting these particular final states. While 
some combinations of cuts in these quantities are more efficient for the signal than the back- 
grounds, there is no overall gain in the signal significance, further, to the cost of a reduced 
signal absolute rate. (Results are here shown for just one signal mass value, yet they are 
typical also for the other masses considered too.) 

4 Conclusions 

The Higgs sector of the NMSSM is phenomenologically richer than the one of the MSSM 
as it has two more (neutral) Higgs states. In this paper, we explored the detectability 
of the lightest CP-odd Higgs state of the NMSSM, a u at the LHC and SLHC, through its 
production in association with a bb pair followed by a\ — » bb decay. We have shown that there 
are some regions of NMSSM parameter space where a rather light a\ state can be discovered 
at very large integrated luminosity using standard reconstruction techniques [2H[25J. This 
is true so long that tan/3 is very large, typically above 30, over an m ai interval extending 
from 20 GeV to 80 GeV, depending on the collider configuration. High 6-tagging efficiency 
is crucial to achieve this, to transverse momenta as low as 15 GeV, ideally, also accompanied 
by di-jet mass resolutions up to 10 GeV. 

Firstly, these results (together with those of Refs. [TTj ITS] ) support the 'no-lose theorem' 
by looking for the (quite possibly resolvable) direct production of a light ai rather than 
looking for its production through the decay hi t 2 —> didi, whose detectability remains uncer- 
tain. Secondly, they enable distinguishing the NMSSM Higgs sector from MSSM one (thus 
contributing to enforcing a 'more-to-gain theorem') since such light CP-odd Higgs states, 
below Mz (which have a large singlet component), are not at all possible in the context of 
the MSSM. 

In reaching these conclusions we have, first, sampled the entire parameter space of the 
NMSSM to extract the regions which are interesting phenomenologically and, then, sampled 
five benchmark points in the NMSSM parameter space (see the captions to Figs. 3-7), all 
taken at large tan/3, for which we have extracted the detector performances and collider lu- 
minosities required for discovery. The encouraging results obtained should help to motivate 
further and more detailed analyses. 
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;ure 1: The rates for o(gg — >• bbai) BR(ai —> bb) as a function of A, k, tan/3, /i eff , A\, A 
ai — > bb) and m ai . 8 
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Figure 2: The BR(ai — > bb) as a function of the CP-odd Higgs mass m ai and of the BR(ai 
77) • 
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Figure 3: The differential cross-section in di-jet invariant mass rrijj, after the cuts in ([T]), for 
signal (with m ai =19.98 GeV) and backgrounds, the former obtained for the benchmark point 
with A = 0.075946278, re = 0.11543578, tan/3 = 51.507125, /i = 377.4387, A x = -579.63592 
and A K = -3.5282881. 
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Figure 4: The differential cross-section in di-jet invariant mass rrijj, after the cuts in ([TJ, for 
signal (with m 0l =35.14 GeV) and backgrounds, the former obtained for the benchmark point 
with A = 0.091741231, re = 0.51503049, tan/3 = 38.09842, /i = 130.56601, A x = -720.88387 
and A K = -5.3589352. 
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Figure 5: The differential cross-section in di-jet invariant mass rrijj, after the cuts in ([TJ, for 
signal (with m ai =46.35 GeV) and backgrounds, the former obtained for the benchmark point 
with A = 0.14088263, re = 0.25219468, tan/3 = 50.558484, /i = 317.07532, A A = -569.60665 
and A K = -8.6099538. 
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Figure 6: The differential cross-section in di-jet invariant mass rrijj, after the cuts in ([TJ, for 
signal (with m 0l =60.51 GeV) and backgrounds, the former obtained for the benchmark point 
with A = 0.17410656, re = 0.47848034, tan/3 = 52.385408, /i = 169.83139, A x = -455.85097 
and A K = -9.0278415. 
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Figure 7: The differential cross-section in di-jet invariant mass rrijj, after the cuts in ([TJ, for 
signal (with m ai =80.91 GeV) and backgrounds, the former obtained for the benchmark point 
with A = 0.10713292, re = 0.13395171, tan/3 = 44.721569, fi = 331.43456, A x = -418.13018 
and A K = -9.7077267. 
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Figure 8: The significances S/y/B (left) and total event rates S (right) of the gg,qq — > 
bbai — > bbbb signal as functions of the integrated luminosity for 10 GeV (top) and 5 GeV 
(bottom) di-jet mass resolutions. 




Figure 9: The differential cross-section in di-jet (pseudorapidity-azimuth) separation, after 
the cuts in (jTJ), for signal (with m ai =46.35 GeV) and backgrounds, the former obtained 
for the benchmark point with A = 0.14088263, k = 0.25219468, tan/3 = 50.558484, /i = 
317.07532, A\ = -569.60665 and A K = -8.6099538. On the left the absolute normalisations 
are given while on the right those to 1. 
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Figure 10: The differential cross-section in jet transverse momentum, after the cuts in ([I]), 
for signal (with m ai =46.35 GeV) and backgrounds, the former obtained for the benchmark 
point with A = 0.14088263, k = 0.25219468, tan/3 = 50.558484, /z = 317.07532, A x = - 
569.60665 and A K = -8.6099538. On the left the absolute normalisations are given while on 
the right those to 1. 
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